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I/") , A condensate initially prepared at finite temperature evolves under external time-dependent 

' perturbation into a time-dependent mixed state. In these notes I use number-conserving time- 

dependent Bogoliubov theory to derive probability distribution for different outcomes of density 
measurement on the time-dependent excited state. 
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I. INTRODUCTION 

Quantum measurements on Bose-condensed systems can give quite unexpected results. For example, in the classic 
paper by Javanainen and Yoo [1] a density measurement on a Fock state | AT/2, N/2) with N particles equally divided 
between two counter-propagating plane waves e ±lx reveals an interference pattern p(x\tp) ~ cos 2 (x — tp) with a 
phase tp chosen randomly in every realization of the experiment. The Fock state has a uniform single particle density 
distribution, but its measurement unexpectedly reveals interference between the two counter-propagating condensates. 
The Fock state is a quantum superposition over Y-particle condensates with different relative phases tp in their wave 

S functions [2], \N/2, N/2) ~ J dip \N : e +i ( x ~ tp > + e~' i ( a: ~ ¥ ')), but every single realization of the experiment reveals such a 
density distribution as if the state before the density measurement were one of the condensates \N : e +^-v) +e -*{*-v)) 
with a randomly chosen phase <p. This effect is best explained [2] when the density measurement, which is a destructive 
measurement of all particle positions at the same time, is replaced by an equivalent sequential measurement of one 
position after another. With an increasing number n of measured positions a quantum state of the remaining N — n 
particles gradually "collapses" from the initial uniform superposition over all phases to a state with a more and more 
localized phase tp. For a large N a measurement of only a small fraction -S* <C 1 of all particles practically collapses 
the state of remaining N — n particles to a condensate with definite phase tp. 

A lesson from this instructive example [1,2] is that quantum measurement on an Y-particle state with highly 
occupied single particle modes "collapses" the state to a definite condensate with a definite condensate wave function 
4>{x). The question is: what is the probability distribution for different measurement outcomes <fi(x)7 As the set of 
condensates is not an orthonormal basis this is not a trivial question. 
' In Ref. [6] we derived this probability distribution in the framework of the time-dependent Bogoliubov theory at 
zero temperature. At zero temperature a condensate initially prepared in its Y-particle ground state evolves under 
external time-dependent perturbation into a time-dependent excited state. The excited state is a time-dependent 
Bogoliubov vacuum i.e. at any time t there exists a complete set of quasiparticle annihilation operators for which the 
excited state is a vacuum. In Ref. [6] it was shown that the time-dependent vacuum has a simple diagonal structure 
which directly leads to a compact gaussian probability distribution for different condensate wave functions 4>{x). As 
the case of zero temperature is covered in Ref. [6], in these notes I describe the general case of finite temperature 
when the initial state is a condensate in equilibrium with a thermal cloud of atoms. I derive gaussian probability 
distribution for <j)(x) at any time t when the external perturbation drives the initial thermal state into an excited 
^ mixed state. 
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II. N-CONSERVING BOGOLIUBOV THEORY 



Number conserving Bogoliubov theory [7] is a quadratic approximation to the second quantized Hamiltonian (in 
trap units) 



H = dx 



-d x &d x $ + ix 2 *^ + V(t, x)&i> + -g&&M 



(1) 
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Here ^(x) is the bosonic annihilation operator, V(t,x) is the external perturbation potential, and g is strength of 
contact interaction between atoms. Here and in the following I will use one-dimensional notation but all equations 
can be generalized by the simple replacement x — > x. The annihilation operator is split into condensate and non- 
condensate part 

*(a;) = a 4>o(x) + Sip{x) . (2) 

It is assumed that most atoms occupy the condensate mode 4>o(x)- Equation (2) is substituted to the Hamiltonian 
(1) and then the Hamiltonian is expanded in powers of the fluctuation operator Sip, see Ref. [7]. 

Many experiments on dilute atomic condensates can be clearly divided in two steps: as a first step a condensate 
is prepared in its ground state and then in the second step an external potential V(t, x) is applied to manipulate 
with the condensate wave function. Generic examples are phase imprinting of dark solitons [5] , atomic intcrferometry 
[3], or generation of shock waves in Bose- Einstein condensates [4]. At finite temperature the initial state before 
the manipulation is a thermal state including thermal excitations above the iV-particle ground state. The initial 
condensate wave function <j> solves the stationary Gross-Pitaevskii equation 

M) = ~dl4>o + 7^x 2 4>o + g\<f>o\ 2 <h • (3) 

In Bogoliubov approximation the ground state is Bogoliubov vacuum |0&) which can be written as a gaussian super- 
position over condensates [8] 
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Here the state \N : (f>) is a condensate of N atoms in the normalized condensate wave function , 2 The Bogoliubov 

\/(4>\<t>) 

modes u m and v m are eigenmodes of the stationary Bogoliubov-de Gennes equations 

u m u m = -\^ 2 x u m + ^x 2 u m + 2g\(f) \ 2 u m + gcj)lv m , 

-u! m v m = ~dlv m + ^x 2 v m + 2g\4> a \ 2 v m + g (<j)* ) 2 u m . (5) 

Numerical solution of these equations gives a finite number of modes M. At finite temperature the initial state is a 
thermal state p(0) with thermal quasiparticle excitations. The thermal state is also a gaussian state [8] 

p(0) = y d 2 b L J d 2 b R e-^-ibnbn+ble-^bn 

1 M \ I 1 M 

* m=l / \ * m=l 

Here b* L b L = ^J>L, m and b\e~^b R = £^=i b^ m e-^b R , m . 

In Bogoliubov theory the initial thermal state evolves under external perturbation V(t, x) into an excited state 
p{t) which has the same form as the initial /5(0) in Eq.(6) but with time-dependent Bogoliubov modes u m (t, x) and 
v m (t, x) which solve time-dependent Bogoliubov-de Gennes equations 

id t u rn = -\d 2 x u m + \x 2 u m + 2g\(j) \ 2 u m + geftv. 
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-id t v m = ~dlv m + ^x 2 v m + 2g\(j) \ 2 v m + g (<f>* a ) 2 u m (7) 

with initial conditions being the eigenmodes of the stationary BdG equations (5). The time-dependent condensate 
wave function <j> (t,x) solves the time-dependent Gross-Pitaevskii equation 

id t (f> = ~dl<j)o + ^x 2 ^ + g\4>o\ 2 (j)o ■ (8) 
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III. PROBABILITY DISTRIBUTION FOR OUTCOMES 



As mentioned before, density measurement is "collapsing" A-particlc state to a Bose-Einstein condensate. The aim 
of the measurement theory is to provide probability distribution for different condensate wave functions <j>. In the 
present context of Bogoliubov theory it is convenient to split possible condensate wave functions into the condensate 
part and the non-condensate part: 4> = (f> + -^8<f>. The aim is to find gaussian probability distribution for 5(f) in the 
gaussian state p(t). Ideally the gaussian distribution would be fully determined by the following equalities between 
second order correlators of the gaussian 8<j>(x) and second order correlators of the field operators: 

6(j)*(x)6(f)(y) = (5fo(x)6ip(y)) = ^n m u* m {x)u m (y) + (1 + n m )v m (x)v* m (y) , (9) 

m 

WW = (S^(x)S^(y)} = £ (x)v m (y) + (1 + n m )u m (x)u* m (y) , (10) 

m 

6<j)(x)8<j)(y) = (6i>(x)6ip(y)) = ^ n m v* m {x)u m {y) + (1 + n m )u m {x)v* m {y) , (11) 

m 

5<j)* (x)S(f>* (y) = (6$ (x)6$ (y)) = ^ n m u* m (x)v m (y) + (1 + n m )v m (x)u* m (y) . (12) 

m 

Here the most right hand sides follow from the Bogoliubov theory [7]. n m — — 1) _1 is average number of 

thermally excited Bogoliubov quasiparticles in the initial state. Unfortunately, because of the non-zero commutator 
[8^j(x),S^(y)] = S(x — y) — (f>l(x)(j)o(y), the first two conditions cannot be satisfied simultaneously. I replace them 
with the condition 

WWW) = (SfixWiv))* = ±(6ft(x)84>(y) + 64,(y)64>\x)) . (13) 
It is convenient to expand the fluctuation as 

oo 

8<j)(x\z) = ^ z a(f>a(x) (14) 
a=l 

in the ortonormal basis of the eigenmodes <j) a of the reduced single particle density matrix 

M 

(S^(x)S^(y)) = SN a <f a (x)<f> a {v) ■ (15) 



a=l 



Here the left hand side is given by Eq.(9). The right hand side is obtained after diagonalization of the hermitean 
operator on the left. The real eigenvalues SN a are average occupation numbers of the corresponding non-condensate 
modes <p a . The correlators (12) after the semiclassical approximation (13) determine the matrix of correlators of the 
complex gaussian random variables z a : 



z a z (3 z a z f3 \ _ ( D aj 3 C a p 
Z^Zp Z a Z* J \ C* a(i D a p 



(16) 



Here the M x M matrices on the right hand side are 



D afj = i (U* m n m Up m + V* m (l + n m )V (im + U am (l + n m )U^ m + V am n m Vp m ) , (17) 
c oc0 = U* am n m Vp m + V* m (l + n m )U* 0m , (18) 

with the matrix elements U am — {4> a \u m ) and V am = {4> a \vl n ). Replac ing Zq-'s with real coordinates, z a — x a -\- iya, 
we get a real symmetric matrix of correlators 



x a X{i x a yf3 \ = 1 / Re D a/3 + Re C Q/3 Im C afs 

y^xp TMJp J 2 V Im c a(i Re D a p - Re C a/3 



(19) 



and a condition that Im D a p = - a good test of the correctness of the calculations. Diagonalization of the correlation 
matrix (19) gives eigenvalues A s > with s = 1, ..,2M. Corresponding eigenvectors are columns of an orthogonal 
matrix O. The eigenvectors define convenient parametrization of the gaussian fluctuation as 



3 



M M 2M 2M 

= ^2 z a <j> a (x) = ^2 4> a {x) ^2 (O a . s + iO M+a ,s) q s = ^2 ®s(x)q s (20) 

a— 1 a— 1 s— 1 s— 1 

with independent real gaussian random variables q s of zero mean and variances g 2 — A s . However, this is not the end 
of the story yet. 

As a result of the semiclassical approximation in Eq.(13) averages like e.g. average density of depletion 5<fi* (x)6<f)(x) 
are divergent because there is infinite number of unoccupied modes (j> a (x), every one of them contributing to this 
depletion density a term i^>* (x)<f> a (x). In stochastic averages like 5<p* (x)S(f>(x) average occupation numbers of modes 
(j> a (x) seem to be 5N a + i instead of the correct 5N a . This artifact of the semiclassical approximation can be corrected 
by introducing to Eq.(20) of regularization factors: 

m / SN \ 1/2 M { SN \ 1/2 2M 2M 

$<P(x)\ Icg = ^2z a i SN fr j <f> a (x) = J2 ( gN 1 l ) Mx)^2(O a ,s +iO M + a ,s)q s = ^2<S> s (x)q s (21) 

a=l ^ a 2 ' Q= l ^ a 2 ' s =l s=l 

/ \ !/ 2 

As expected in semiclassical approximation, the regularizing factors ( g ^_i ) are approximately 1 for the highly 
occupied modes with SN a 3> 1 which dominate in the density distribution, but at the same time they remove the 
divergence coming from the infinity of unoccupied modes. The wave functions $ s (x) are in general neither normalized 
nor orthogonal, except in the quantum limit of zero temperature, see the proof in Ref. [6]. 

IV. CONCLUSION 

In conclusion, a recipe to simulate density measurement on the time-dependent excited thermal state has the 
following steps: 

• Solve stationary Gross-Pitaevskii and Bogoliubov-de Gennes equations (3,5) to provide initial conditions for 
4>o(t,x), u m (t,x) and v m (t,x), and the initial quasiparticle frequencies ui m . 

• Solve time-dependent Gross-Pitaevskii and Bogoliubov-de Gennes equations (8,7) with respect to (j)o(t,x), 
u m (t,x), and v m (t,x). 

• Diagonalize the reduced single particle matrix (15) to get its non-condensate eigenmodes <fi a with their average 
occupation numbers SN a . 

• Build the matrices D a p and C a p in Eqs.(17,18), and then the real symmetric correlation matrix in Eq.(19). 

• Diagonalize the correlation matrix in Eq.(19) to get its real eigenvalues \ s and correposnding eigenvectors O as . 

• Build the regularized modes $ s according to their definition implicit in Eq.(21): 

™ / X AT \V2 



a=l V" iV " + 2 / 



(22) 



• Choose independent real random variables g s 's from their gaussian distributions of zero mean and variance 
(j 2 = A s , and then combine the chosen q's into condensate density 



p(x\q) 



2M 



N<f> (x) + J2q s <5> s (x) 



(23) 



The p(x\q) defines a family of all possible density measurement outcomes with a gaussian probability distribution for 
different q's. 

In the limit of zero temperature this general recipe coincides with the recipe derived by different methods in Ref. 
[6]. At zero temperature the wave functions become <& a (x) ~ (/>a(z) for a = 1, M and zero otherwise (here the ~ 
means equality up to a phase factor). Corresponding variances are g 2 = SN a . 
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